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Abstract. We report on the measurement of D meson production from the analysis 
of their hadronic (D° — > Kir) and semileptonic (D —> fi + X, D ^ e + X) decays in 
y/sNN=200 GeV Au+Au collisions. The transverse momentum (pt) spectra and the 
nuclear modification factors for D° and for electron/muonj from charm semileptonic 
decays will be presented. The differential cross section da/dy is found to be consistent 
with the number of binary scaling. The blast-wave fit suggests that the charm hadron 
freeze out earlier than other light flavor hadrons. 

PACS numbers: 25.75.Dw, 13.20.Fc, 13.25.Ft, 24.85. +p 
1. Introduction 

In relativistic heavy-ion collisions, charm quarks are predicted to lose less energy 
compared to light quarks in the partonic matter due to the "dead-cone" effect [UI21IH]- 
The pt distributions and nuclear modification factors of D mesons and of single electrons 
from charmed hadron decay will be vital to study the physics in relativistic heavy-ion 
collisions 4J. Charm quarks are believed to be produced at early stages via initial gluon 
fusion and their production cross section can be evaluated by perturbative QCD [Oj. 
Study of the binary collision (Nun) scaling properties for the charm total cross section 
among p+p, d+Au to Au+Au collisions can test if heavy-flavor quarks are produced 
exclusively at initial impact [TP] . Due to the heavy mass of charm quark, charmed 
hadrons might freeze out earlier than light flavor hadrons. Their flow velocity might be 
a good indicator of light flavor thermalization occurring at the partonic level [HJ El HI IHj ■ 



2. Analysis and Results 

The data used for this analysis methods were taken with the STAR experiment during 
the y/sj^=200 GeV Au+Au run in 2004. A total of 13.3 and 7.8 million 0-80% minimum 

f For the full list of STAR authors and acknowledgements, see appendix 'Collaboration' of this volume. 
X The word "muon" refers to "electron" refers to electron/positron throughout these 

proceedings. 
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Figure 1. Left panel: Upper part: Particle mass squared distribution (m 2 = 
(p/P/j) 2 ) from the TOF after TPC dE/dx selections. A clear muon mass peak is 
observed and the primary pion candidates are shown as the right peak. Bottom part: 
Primary particle DCA (red line) and muon DCA distributions from pion kaon weak 
decayed background after TPC dE/dx and TOF m 2 selections from HIJING simulation 
through realistic STAR detector configuration (blue line). Right panel: The pT spectra 
of D° and electron/muon from semileptonic decays in Au+Au and d+Au collisions. 
Dashed curves are from the power-law fits by combining above three measurements, 
while solid curves are from the blast-wave fits by combining D° and muon spectra with 
non-photonic electron spectra at pr < 2 GeV/c. 

bias Au+Au events were used for the D° reconstruction and the Time-of-Flight (TOF) 
single muon/electron analysis, respectively. About 15 million top 12% central Au+Au 
collision events were also used for TOF single muon and electron analysis. 

The D° mesons (px < 3 GeV/c) were reconstructed through their decay D° — ► 
K~-k + (D° — > K + tt~) with a branching ratio of 3.83%. The D° yields were obtained 
from the invariant mass distributions of kaon-pion pairs after mixed-event background 
subtraction. Analysis details can be found in Ref. [TTj . 

Inclusive electrons can be identified up to pr = 5 GeV/c by using a combination 
of velocity (/?) from the TOF and ionization energy loss (dE/dx) measured in the 
STAR Time Projection Chamber (TPC). To measure the photonic electron spectra, the 
invariant mass and opening angle of the e + e~ pairs were constructed from an electron 
(positron) in the TOF at lower px (< 1.2 GeV/c) or in the TPC at higher p? (1.2 < 
Pt < 5 GeV/c) combined with every other positron (electron) candidate reconstructed 
in the TPC OUB]. The non-photonic single electron spectra (0.9 < px < 5 GeV/c) 
from charm semileptonic decays can be extracted from the inclusive electron spectra 
subtracted by photonic background. Detailed analysis can be found in Ref. [TTJ[T3|. 

However, due to the large combinatorial background in charmed hadronic decay 
channels and the overwhelming photon conversions in the detector material, we 
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Figure 2. Left panel: Charm hadron freeze-out temperature Tf a versus maximum 
now velocity f3 m . Right panel: Charm cross section at mid-rapidity compared with 
pQCD calculations in d+Au and Au+Au collisions. 



conducted the measurement of single muon spectra at low pt (0.17 < px < 0.25 GeV/c) 
in both 0-80% minimum bias and top 12% central Au+Au collisions at y/s NN = 200 
GeV. Muon identification was provided by dE/dx from the TPC and velocity from 
the TOF. Background muons from pion/kaon weak decays were subtracted using the 
distribution of the distance-of-closest-approach (DCA) to the collision vertex. The single 
muon raw yield was obtained from the fit to muon DCA distributions by combining the 
background DCA distributions and the primary particle DCA distributions [TO} I14j. 

The spectra of D° and muons/electrons from charmed hadron decays are shown in 
the right panel of Fig. ^ The electron yield relies more on D° spectrum shape with 
a wide range of power-law parameters n and (pt) at fixed charm total yield. A factor 
of 8 variation of electron yield integrated above pt of 1.0 GeV/c was seen, while muon 
yield only changes within ±15% jTHj. We therefore conducted a measurement of muons 
at low pt in order to constrain the charm cross section. The (pt) and n can be derived 
from a power-law fit to the D° px distributions and decayed lepton spectra shape, while 
the freeze-out temperature Tf Q and flow velocity Pt can be derived from a blast-wave fit 
to the D° p T distributions and decayed lepton spectra shape (px < 2 GeV/c), see the 
right panel of Fig. The la contour from the combined blast- wave fit with a quadratic 
sum of stat. and sys. errors for the spectra shows the lower limit of charm hadron Tf Q 
(> 140 MeV) and small but non-zero p T (0.21 ± 0M(stat.) ± 0.07 (sys.)). This result 
suggests that charm hadron seems to freeze out earlier than other hadrons, see the left 
panel of Fig- El The charm cross sections at mid-rapidity (da/dy) can be obtained from 
the average of the two fits by combining the three measurements covering ~ 90% of 
the kinematics. They are 301 ±44(stat.) ± 67(sys.)fj,b in 200 GeV d+Au collisions jllj . 
267 ± 19 ± 49/z& in 200 GeV 0-80% minimum bias Au+Au collisions and 283 ± 12 ± 39/jb 
in 200 GeV top 12% central Au+Au collisions. Within error bars, the measured charm 
differential cross sections are found to be consistent with the number of binary scaling, 
see the right panel of Fig. El In addition, the measured cross sections are larger than the 
pQCD prediction by a factor of 5 [9 ! Note that the systematic errors are the dominant 
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ones. In the fitting procedure, the statistical and systematic errors were summed up 
quadratically. 



3. Conclusions 

The measurement of charm production from analysis of D° — * Ktc, muon, and electron 
channels in both minimum bias and central Au+Au collisions at RHIC was reported. 

The transverse momentum spectra from non-photonic electrons are strongly 
suppressed at 1 < px < 5 GeV/c in Au+Au collisions relative to that in p+p and 
d+Au collisions [4 . For electron with px > 2 GeV/c, corresponding to charm hadron 
with pt > 4 GeV/c, the suppression is similar to those of light baryon and meson 
hadrons 

Detailed model-dependent analysis of the electron spectra with px < 2 GeV/c 
suggests that charm hadrons have a different freeze-out pattern from the more copiously 
produced light hadrons. The blast-wave fits show that charm hadrons seem to freeze 
out early, and have similar freeze-out temperature as multistrangeness hadrons (<p, Q) 
but with smaller collective velocity. 

Charm cross sections at mid-rapidity are extracted from a combination of the three 
measurements covering ~ 90% of the total yield at mid-rapidity. The cross section is 
found to follow binary scaling, which is a signature of charm production exclusively at 
the initial impact. This supports the assumption that hard processes scale with binary 
interactions among initial nucleons and charm quarks can be used as a probe sensitive 
to the early dynamical stage of the system. 
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